Excitations from a chiral magnetized state of a frustrated quantum spin liquid 
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We study excitations in weakly interacting pairs of quantum spin ladders coupled through ge- 
ometrically frustrated bonds. The ground state is a disordered spin liquid, that at high fields is 
replaced by an ordered chiral helimagnetic phase. The spectra observed by high- field inelastic neu- 
tron scattering experiments on the prototype compound SUI-CU2CI4 are qualitatively different from 
those in the previously studied frustration- free spin liquids. Beyond the critical field He = 3.7 T, 
the soft mode that drives the quantum phase transition spawns two separate excitations: a gap less 
Goldstone mode and a massive magnon. Additional massive quasiparticles are clearly visible below 
ifc, but are destroyed in the ordered phase. In their place one observes a sharply bound excitation 
continuum. 
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I. INTRODUCTION 

Gapped quantum-disordered antiferromagnets (AFs), 
also known as "spin liquids" , have become prototype ma- 
terials for the study of Bose-Einstein condensation and 
related phenomena i^^i^i^i^i^ By virtue of the Zeeman ef- 
fect, the magnetic field directly affects the chemical po- 
tential for the relevant lowest-energy triplet excitations. 
The latter excitations can be viewed as bosons with hard- 
core repulsion. Exotic quantum critical points are real- 
ized when the gap energy for one member of the triplet 
is driven to zero at some critical field i^c, and a macro- 
scopic number of those magnons get incorporated into the 
ground state. At H > He one typically observes a mag- 
netically ordered state. However, due to the strong inter- 
actions the corresponding excitation spectrum is nothing 
like that of conventional ordered AFs.^ 

Usually dit H > He the magnon branch that goes soft 
at the quantum critical point is replaced by a gapless col- 
lective Goldstone modeA'^^ For all models and prototyp- 
ical materials studied to date, the two other members of 
the original excitation triplet survive as sharp gap excita- 
tions. Thus, the low-energy spectrum remains dominated 
by three long-lived quasiparticles. In the present work we 
study the spin dynamics in a gapped quantum AF with 
strong geometric frustration of magnetic interactions. In 
a stark contrast to the behavior of frustration-free spin 
liquids, we find that above He the soft mode gives rise to 
two distinct excitation branches, while the higher-energy 
gapped magnons become unstable and are replaced by a 
broad excitation continuum. 

Our prototype material, Sul-Cu2Cl4,^^^^^ is a quasi- 
one-dimenisonal Heisenberg AF with a singlet ground 
state and a spin gap of Aq = 0.52 meV. The corre- 
sponding network of S = 1/2 Cu^+ spins is best de- 
scribed as an array of 4-leg spin tubes with dominant 
AF nearest-neighbor coupling along the legs and sev- 
eral weaker rung interactions of comparable strength. 
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FIG. 1: A schematic representation of the magnetic interac- 
tion topology in SUI-CU2CI4. Ji denotes the "leg" coupUng, 
and J2, J3 and J4 are the three distinct "rung" interactions. 



The tube's legs run along the crystallographic c axis of 
the triclinic PI structure. The corresponding interac- 
tion topology is illustrated in Fig. [TJ Each spin tube 
consists of two spin ladders, Ji and J3 being the leg 
and the rung exchange constants, respectively. The two 
ladders are coupled via exchange constants J2 and J4, 
that are in obvious geometric frustration with Ji. Even 
though the ground state is disordered, frustration en- 
sures that dynamic one-dimensional spin correlations in 
SUI-CU2CI4 are peaked at incommensurate wave vectors 
lo = 0.5 - 5, 1.5 + J, J = 0.022(2)^ The dispersion of the 
gap excitation along the legs of the tubes can be approx- 
imated as 



(1) 



where vq ~ 14 meV is the spin wave velocity. Any disper- 
sion transverse to the c axis is undetectably small. Pre- 
viously we have demonstarated that the application of a 
magnetic field exceeding He = 3.7 T destroys the spin 
liquid and induces an ordered helimagnetic state with an 
incommensurate propagation vector.^ The incommensu- 
rability along the c axis exactly corresponds to the min- 
imum of the dispersion at zero field. The present work 
deals with the spin dynamics on either side of this field- 
induced phase transition. 
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II. EXPERIMENTAL PROCEDURES 

Inelastic neutron scattering measurements were per- 
formed on the same fully deuterated single crystal sam- 
ples as those used in Ref. 8. Two series of experiments 
were carried out at the V2-Flex and TASP cold neutron 
spectrometers at HMI and PSI, respectively. In Setup 1 
we used Ef = 3.7 meV fixed final energy neutrons, while 
the sample environment was a 14.5 T split-coil cryomag- 
net with a dilution refrigerator insert. Data in Setup 
2 were collected with Ef = 3 meV and a similar 9 T 
setup. In both cases we utilized pyrolitic graphite (PG) 
vertically focused monochromators and horizontally fo- 
cused analyzers, as well as Be filters positioned after the 
sample. The field was in all cases applied along the b 
axis and scattering data were collected in the (/i, 0, /) 
reciprocal-space plane. Most of the data were measured 
in constant-g scans at the wave vector transfers where 
the previously measured ID magnon dispersion is a min- 
imum: g=(/i, 0, 0.48). Momentum transfers along the a* 
axis were chosen to optimize wave vector resolution along 
the c* direction. The background was measured at the 
wave vectors (/i, 0,0.6), where no magnetic contribution 
is expected, due to a very steep dispersion along the crys- 
tallographic c* direction. These background scans were 
fit to a constant plus a Gaussian function centered at 
zero energy transfer, to account for fast neutron back- 
ground and elastic incoherent scattering from the sample 
and sample environment, respectively. The resulting fit- 
ted background function was subtracted from the signal 
scans. 
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III. EXPERIMENTAL RESULTS 

Typical background-subtracted scans collected at < 
He are shown in Fig. [2l As the external magnetic field is 
increased, the single peak seen at = (Ref. 1^, Fig. 3) 
splits into three components. The gap energy of the cen- 
tral and most intense component (mode 2) is field inde- 
pendent. The energies of modes 1(3) decrease (increase) 
with the field, as expected for Zeeman splitting of the 

= 1 excitation triplet. The gap in the lower branch 
(mode 1) approaches zero energy diS H ^ He. 

For magnetic fields exceeding He = 3.7 T (Fig. [3]) 
the measured excitation spectrum undergoes some dras- 
tic changes. As can be seen in Fig. [3^, just above the 
critical field, at = 4 T, mode 3 becomes very weak 
and is, in fact, barely visible. It is totally absent for all 
other values of magnetic field applied in our experiments 
(5 T, 8 T, 9.5 T and 13.5 T). The gap energy for mode 
2 starts to increase with increasing H. Finally, mode 1, 
being critical H = Hc^ re- acquires an energy gap gap 
at higher fields (Fig. [3]b,c). This gap, however, increases 
rather slowly with the field not exceeding 0.3 meV even 
at the highest attainable field of 13.5 T. 



FIG. 2: Field dependence of background-subtracted inelastic 
scattering at the ID AF zone-center (0.5, 0, 0.48), collected 
at H < He and T ^ 70 mK. Solid lines are model cross sec- 
tion fits, as described in Section |Vl Shaded areas are partial 
contributions of three separate excitation branches. Inserts: 
Dispersions of the three field-driven excitation branches in 
SUI-CU2CI4 and evolution of the calculated FWHM resolu- 
tion ellipsoids plotted in projection onto the (/, hu) plane. 



IV. THEORETICAL CONSIDERATIONS 

The re-opening of the spin gap dX H > He is, at 
the first glance, puzzling. Indeed, in those spin liquids 
that undergo similar field-induced ordering transitions 
the gap re-appears in the high field phase due to mag- 
netic anisotropy effects. In that case the transition is of 
the Ising, rather than the BEG, universality class. How- 
ever, in the case of SUI-CU2CI4, magnetic anistropy is 
negligibly small, as discussed in Ref. 9. Moreover, since 
the helimagnetic ordering occurs at an incommensurate 
wave vector, one can rigorously prove that there must 
be a gapless "sliding mode" ("phason"). This remains 
true even if anisotropy is present. The sinusoidal spin 
spiral that is the ground state in the isotropic case, will 
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FIG. 3: Background-subtracted neutron spectra measured in 
SUI-CU2CI4 at the ID AF zone-center (0.5, 0, 0.48) in vari- 
ous applied magnetic fields above the critical field, H > He, 
and 150 mK. The lines, shaded areas, and graphical repre- 
sentations of the resolution functions are as in Fig. (2] The 
shaded area in the inserts indicates a continuum scattering, 
with the corresponding solid line depicting the lower bound 
of the continua. 



be distorted by the anisotropy to form a soliton lattice. 
That soliton array, however, being incommensurate, will 
retain the ability to slide freely relative to the underlying 
crystal structure. 



To fully understand the nature of the spin gap in Sul- 
CU2CI4, we shall develop a field-theoretical description 
of the low-energy spin correlations in a pair of magne- 
tized spin ladders coupled with geometrically frustrated 
interactions. 



A. Inter-ladder interactions 

We start with writing down the spin Hamiltonian for 
the frustrated spin tube geometry shown in Fig. [1] 

n;a=l,2 

+ ^ [J4Sl^Sl + J2S^^iS^+ 

n 

+ J2TI T2 + (2) 

Here a = 1,2 labels the two structural ladders, while 
Sn an are spin operators on the corresponding legs. 
We direct the reader's attention to the fact that the ratio 
of spin gap Aq to the magnetic excitations bandwidth vq 
in SUI-CU2CI4 is much smaller than in most known spin 
ladders. This is likely to be a consequence of geometric 
frustration. Whatever the cause, it makes this particular 
system ideal for application of the field theory methods 
developed in Ref. 10. 

We shall treat Hamiltonian ([2|) as one describing two 
weakly interacting ladders. Assuming Ji ^2,-^3,-^4, 
there are two cases to consider: (i) J3 ^ | J2 — J4I = J±. 
Here one ladder is composed of spins 81^2 and the other of 
spins Ti^2- Such ladders have spectral gaps A ~ J3. The 
second sum in ([2]) represents a frustrated interaction and 
in the gapped phase can be considered as a perturbation. 
This is justified by the fact that the term with J3 has 
smaller scaling dimension than J2 , J4 hence generating 
a stronger coupling, (ii) Alternatively, one may consider 
I J2 — J4 1 ^ J3 = J±. Then the two ladders are composed 
of S and T spins, respectively. The spin gap is of the 
order of IJ2 — J4I. In either scenario, the inter-ladder 
interactions J± are irrelevant at small magnetic fields, 
yet become relevant when the external field drives the 
single-ladder gap to zero. 

To obtain the continuum limit of the lattice model ([2]), 
we approximate the spin operators as follows: 

^M^(x) + (-l)^n^(x), a = 1,2 
T^^L'^(x) + (-l)^F(x), a = 1,2 (3) 

where 1, n and L, M are operators with scaling dimen- 
sions 1/2 and 1 respectively. It turns out that only 
antibonding combinations of the staggered magnetiza- 
tions become critical in high magnetic fields. There- 
fore only such combinations will generate relevant inter- 
ladder coupling dit H > He. They are defined as = 
- P (a = 1,2) and Ni = ni - n2,N2 = li - I2, for 
cases (i) and (ii), respectively. 

Substituting these definitions into ([2]) and keeping only 
the antibonding staggered magnetization terms, we ob- 
tain the following inter-ladder interaction: 

V = jjdx [Ni5,N2 - N2a,Ni] , (4) 

where 7 ~ J_l/Ji. 
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B. Ginsburg-Landau model 

Interaction (|4]) are easier to treat in the spin-liquid 
phase when the spectrum is gapped. Here we shall follow 
Affleck-^ -"^ who suggested to describe the low lying modes 
of spin ladder by an effective Ginzburg-Landau model 
with the Lagrangian density 



2\2 



(5) 



where $ ~ N, and the bare parameters i), A are of the 
same order as the velocity and gap in the actual single- 
ladder spectrum. It is straightforward to show that the 
introduction of the interaction ^ generates two modes 
with the spectrum 



Ei 



(6) 



where are the renormalized gap energy and spin 
wave velocity, respectively, and 5 ^ j. The incommen- 
surate gapped dispersion is fully in agreement with the 
existing data on SUI-CU2CI4. 

When the magnetic field exceeds the spin gap the 
transverse components of the staggered magnetizations 
become critical. Interaction @ becomes relevant, but its 
form is highly non-standard from the field theory point 
of view since the corresponding interaction density has 
nonzero conformal spin (that is its left and right con- 
formal dimensions do not coincide). Attempts to treat 
such interactions using standard field theory methods 
have met with mixed success^^. Below we shall adopt 
a different approach. Again, following Affleck^^, close to 
the critical field we will describe each spin ladder using 
the Landau- Ginzburg theory with the interaction. 
When applying this formalism to our case, one replaces 
= (1 — %lj^ij;)^N^ = ('0,7/^+), where are bosonic 
fields and writes the Lagrangian as 



where cr^ is the Pauli matrix and Q 
interaction density becomes 



where ji ^ H — 
to be large. 



zm 

-^2+5x^i), (7) 
i^c, ^ ^ A, and the coupling g is assumed 



7. Then the 
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[(|r^ip + |r^2p)' + 2|r^ip|r^2n, 



(10) 



since at large Q the oscillatory terms are wiped out. when 
one integrates over x. In terms of rja^ the ground state 
of the resulting model is a ferromagnet. In other words, 
only one species of bosons condenses, e. 771. The 
Hamiltonian for this field is just the bosonic one with a 
point-like repulsion. So the excitation spectrum at small 
wave vectors is linear in wave vector corresponding to 
the phase fiuctuations of the condensate. Bosons of the 
second fiavor remain massive with the energy gap A2 ~ 

Now we have to go back to the original variables. It 
turns out that the magnetic structure observed experi- 
mentally in SUI-CU2CI4 corresponds to the case denoted 
above as (i). For that case we get the following picture 
of the soft spin modes: 



no 



-1? 



N± = (7^lel^^7^^e- 



N 



i(^ie 



\Qx 



(11) 

(12) 



All spin fiuctuations occur in the plane perpendicular to 
the magnetic field. Along the J3 bonds the spins are 
antiparallel. Spins along J2 or J4 are perpendicular to 
each other. The propagation vector along the leg axis is 
27r/o = TT — Q, so /o is close to 1/2. Similarly, for (ii): 



i{r]ie 



iQx 



-iQx\ 



(13) 
(14) 



here the spins are parallel along J2, and perpendicular 
on the J3 bond. 

The model thus qualitatively reproduces the diffraction 
result for the elastic modes, but also predicts that at 
H > He the soft mode splits into two components. One, 
as required, is gapless. The other mode has a gap that 
scales as the square of the transverse ordered moment. 
Below we shall show that this prediction is indeed fully 
consistent with the inelastic data. 



C. Incommensurate order and spectrum at H > He 



We now do the substitution 

^ = l={I^ia,)i^ (8) 

and 

Vii,2 = e±i^^7?i,2, (9) 



V. DATA ANALYSIS 

A quantitative analysis of the inelastic neutron data 
involves fitting the measured scans to a parameterized 
model cross section function, numerically convoluted 
with the 4-dimensional resolution function of the spec- 
trometer. We start with the simpler regime H < He- 
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A. Low fields 



For the vicinity of the minima of the one-dimensional 
dispersion where all scans were collected, the cross sec- 
tion was written in the single-mode approximation as- 
suming three separate excitation branches corresponding 
to Sz = 0, ±1. Since the Zeeman term commutes with the 
Heisenberg Hamiltonian, the spectrum at any H < He 
will be exactly as at = 0, except for a constant shift 
in the energies of the Sz = I and 5'^ = — 1 magnons. 
With this in mind, the dispersion relations for mode a 
{a = 1,2,3) was written as: 



huJa 



+ Ac - Aq. 



(15) 



The spin wave velocity and zero-field gap were fixed at 
vq = 14 meV and Aq = 0.52 meV, respectively, as previ- 
ously measured in zero field. The gap energies Aq, for 
each mode were assumed to be field-dependent. The 
single-mode cross section was then written as: 



dVtdE' 



(X 



^Q!,q)- 



(16) 



Here are are separate intensity prefactors for each of 
the three modes. This model, when convoluted with the 
known instrument resolution, fits all the scans collected 
dit H < He rather well. Such fits are shown in heavy 
solid lines in Fig. O Contributions due to the three com- 
ponents of the spectrum are represented by shaded ar- 
eas. The asymmetric peak shapes with extended "tails" 
on the high-energy side are entirely due to resolution ef- 
fects. 



Ag + vg(|qc-27r/o| + /^)^ 
AolAil Ai^ 
Ao. 



(17) 
(18) 



Typical dispersion relations calculated using this equa- 
tion are the lowest-energy curves plotted in the insets in 
Fig. El 

We shall label the theoretically predicted gapped de- 
scendent of the soft mode with the index a = 4. It 
appears reasonable to approximate its dispersion with 
Eqs. [U and [151 as at < 0. The gap A4 can be accu- 
rately determined from our scans that show a peak at 
the corresponding energy. In contrast, the parameter Ai 
actually represents the velocity at qc ^ Iq. It can not be 
accurately measured due to limitations imposed by the 
experimental resolution. On the other hand, it can not 
be totally disregarded either, as it affects the observable 
high-energy resolution "tail" of the gapless mode. In or- 
der to avoid dealing with an over-parameterized model, 
a compromise was reached by somewhat arbitrarily pos- 
tulating A4 = |Ai|. Indeed, both parameters are zero at 
H = He and are expected to increase in absolute values 
at higher fields. As represented by the heavy solid lines 
in the lower-energy range of the plots shown in Fig. [31 
this model, though somewhat artificial, reproduces the 
experimental data remarkably well. In the same figure 
the light grey shaded areas are partial contributions of 
the gapless and gapped components (modes 1 and 4), 
respectively. 



C. H > He', high energies 



B. H > He', low-energy modes 

Finding an appropriate dispersion relation for the two 
theoretically predicted descendants of the soft mode in 
the regime H > He is not straightforward. The guidance 
provided by the calculations in the previous section is 
limited, as they are only applicable in the direct proxim- 
ity of the ordering vector. 

For the gapless mode it is clear that dit H = He one 
should still be able to employ Eqs. [T5l and [Tl with Ai = 0. 
However, for H > He the dispersion near the ordering 
vector should be linear, with the slope progressively in- 
creasing with H — He. In the strongly ID case like that of 
SUI-CU2CI4, one can expect the dispersion relation to fol- 
low the lower bound of continuum excitations in isolated 
ladders. The latter were investigated in Ref. and can 
be represented by the same parabolic curve as at = 0, 
but offset both in energy and momentum. The "gap" 
Ai becomes negative. The corresponding dispersion and 
neutron cross section are still given by Eqs. [15] and [161 
respectively, but the definition of ficjq in ([T5|) is changed 
to: 



Upon crossing the critical field, the central component 
of the excitation triplet undergoes some drastic changes. 
This mode is the strongest and sharpest feature of the 
spectrum dit H < He (Fig.[2]3), but is replaced by a much 
broader peak above the transition (Fig.[3K), measured us- 
ing in the same experimental configuration). The single- 
mode approximation ([T6|) gives excellent fits to the data 
at low fields, but totally fails to describe the shape of 
the middle mode dit H > He (Fig. [3l dashed curves). 
This leads us to the conclusion that the quasiparticle 
description of this part of the spectrum breaks down 
in the high-field phase. Instead, the excitations are a 
broad continuum of states. By analogy with the Miiller 
ansatz that describes continuum excitations in gapless 
spin chains, we used the following empirical expression 
to approximate the corresponding contribution {a = 2) 
of this spectral feature to the neutron cross section: 

.^^A2\f{q)\\u,-u,2,^r^eiu,-u,2,^). (19) 

In this formula fia;2,q is still given by Eq.[T5l and describes 
the lower bound of the continuum. A new parameteter 
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FIG. 4: Field dependence of the excitation energies at the 
ID AF zone center in SUI-CU2CI4. The points are the data 
from the INS experiments and the sohd hnes are guides for 
the eye. At if > the central excitation branch is replaced 
by a sharp continuum represented by the shaded band. The 
field dependence of the new gapped excitation (above 3.7 T) 
is compared with that of the square of the transverse ordered 
moment, from Ref . 0. 

C characterizes how steeply the intensity falls off with 
energy. With the value ( = 0.30(6) refined at if = 9.5 T, 
this formula yields much better fits to the data at all 
fields dX H > He (Fig. [3l solid curves) than the single- 
mode approximation (dashed curves). 

The highest-energy component of the original triplet of 
gap excitations is also strongly affected by the transition. 
As clearly seen in Fig. [3l its intensity drops dramatically 
and it is no longer observed at if > 4 T. To describe the 
mode at that field however, in our analysis we employed 
an expression similar to that for the middle mode, the 
parameter ^ fixed at the same value. 



VI. DISCUSSION 

Despite fact that our data analysis is based on sev- 
eral somewhat arbitrary assumptions regarding the form 
of the dynamic structure factor, the extracted field de- 
pendencies of the gap energies should be quite reliable. 
Indeed, the gaps are robustly identified by a sharp jump 
in the neutron intensity that occurs in the energy window 
defined by the resolution function. The fact that the cen- 
tral excitation branch is replaced by a sharp continuum 
at if > if c is also a solid experimental finding, as is the 
disappearance of the higher-energy mode. 

The linear Zeeman splitting of the triplet at H < He 
occurs exactly as expected: the system is not ordered 
and the Zeeman term commutes with the Heisenberg 
Hamiltonian. For H > He, the central prediction of 
our Ginsburg-Landau analysis is that the energy of the 
new gapped excitation scales as the square of the trans- 
verse ordered moment. This is indeed consistent with 
the experiment. In Fig. [J] the heavy solid line, and 
the corresponding right-hand-side Y axis, represent the 
field dependence of the latter, as previously measured 
experiment ally To within a constant scaling factor, this 
curve exactly reproduces the field dependence of the re- 
opening gap A4 (if ) measured in this work. 

The dramatic broadening of the higher-energy excita- 
tions H > He is not captured by our hydrodynamic 
model. While a further theoretical study is clearly re- 
quired, we note that the phenomenon resembles the one 
previously observed in the bond-alternating S = 1 chain 
system NTENP.^^ In the latter system it was attributed 
to new magnon decay channels that are opened by the 
lifting of rotational symmetry at if c. A similar mecha- 
nism may at work in SUI-CU2CI4. 



VII. CONCLUSION 



D. Fits to experimental data 

The model cross sections described above were numer- 
ically convoluted withe the spectrometer resolution func- 
tion, calculated in the Popovici approximation^^ They 
were then used in the least-squares analysis of constant- 
Q scans collected at each value of the applied field. The 
resulting fits are shown as the solid lines in Figs. [2] and 
[3l Partial contributions of each spectral component are 
represented by the shaded areas. The insets show the 
dispersion relations of all components calculated using 
the parameters determined in the fitting process. The 
shaded area signifies continuum scattering, with the cor- 
responding solid line indicating the lower bound. In the 
same plot we show the evolution of the experimental res- 
olution ellipse in the course of each constant-Q scan. The 
gap energies obtained in these fits are plotted vs. applied 
magnetic field on Fig. [H (symbols). 



In summary, we suggest that the geometric frustra- 
tion and the coupled-ladder topology of SUI-CU2CI4 are 
responsible for a complex and very unusual excitation 
spectrum in the magnetized state. At the same time, 
triplet low-energy excitations dX H < He are rather typ- 
ical for an isotropic spin liquid. 
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